Spontaneous parity symmetry breaking of 7r-bonds in 2D: 
describing topological insulators in real space 
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The existence of 7r-bonds on topological insulator surfaces is found to be closely related to the 
phenomena of surface conduction and surface band spin polarization. A 7r-bond trimer or 7r-bond 
dimer on the surface can form open conjugated systems that are responsible for the unique surface 
conduction mechanism of a topological insulator. Parity operation in 2D is identified for a 7r-bond 
trimer within a six-fold symmetry coordination and a 7r-bond dimer within a four-fold symmetry 
coordination. Spontaneous 2D parity symmetry breaking is found to be closely related to the 
theoretically predicted ty Berry's phase and the observed surface band spin polarization. The role 
of 7r-bonds on a cleaved Bi2Se3 surface is compared to that for graphene with a 2D band structure 
containing Dirac cones. Similar 7r-bond dimers within a four-fold 2D symmetry coordination can 
also be identified in strained a-Sn of diamond structure as a theoretically predicted topological 
insulator. 

PACS numbers: 73.20.-r; 73.25.+i; 31.10.+z; 33.15.Fm 



I. 



INTRODUCTION 



Topological insulators belong to a novel class of ma- 
terials that exhibits band inversion characteristics due 
to strong spin-orbit interactions, which correspond to 
a unique state of matter for a band insulator with 
surface-only conduction and a spin-polarized surface 
state. Bi2Se3 is a second-generation topological insula- 
tor with a nearly idealized single Dirac cone at the f- 
point. 1 The crystal structure of Bi2Se3 can be described 
in terms of close packed Bi-Se quintuple layers with a 
space group of R3m and van der Waals gaps between the 
Sel-layers^ Topological insulator have long been mod- 
eled using analogous descriptions of topological invari- 
ants, Berry's phase, and the quantum spin Hall effect i£ 
In particular, the defining property of a topological in- 
sulator with a predicted chiral edge mode, i.e., the heli- 
cal Dirac fermions of Z2 topological-order has been con- 
firmed through spin- and angle-resolved photoemission 
spectroscopy (spin-ARPES)4 The success of many so- 
phisticated theoretical predictions and subsequent exper- 
imental verifications since 2006 has led to rapid develop- 
ment in the research conducted on this exotic phase of 
electrons in solids i& 

Electricity conduction by topological insulator surfaces 
is mainly attributed to the band inversion phenomenon 
that has been proposed for heavy elements with strong 
spin-orbit coupling. 6 Band inversion for the partially 
filled p orbitals in Bi and Se may exist under crystal 
field splitting and strong spin-orbit coupling. The crys- 
tal structure of Bi2Se3 can best be described in terms of 
a Bi-Se quintuple- layer unit, where the quintuple layers 
are close-packed with a R3m space group symmetry in 
either a rhombohedral or hexagonal setting*^ Because 
Bi2Se3 has a band gap size of ~0.3 eV in the semiconduc- 



tor range and both Bi/Se elements are close-packed with 
octahedral coordination, we submit that Bi2Se3 may be 
adequately described by a fully hybridized s-p-d orbital 
model, similar to sp 3 hybridization for Si with 3s 2 p 2 , 
rather than by only partially filled p orbitals within an 
ionic model in a point charge crystal field. Recently, we 
proposed that the existence of 7r-bonds and Dirac cones 
in both graphene and the topological insulator surface 
of Bi2Se3 could be the key to explain the unique surface 
conduction mechanism in terms of a local 7r-bond en- 
ergy exchange for graphene and topological insulators ji 
In this study, we report 2D parity symmetry for 7r-bonds 
in topological insulators: the spontaneous 2D parity sym- 
metry breaking may explain experimental observations of 
a spin-polarized gapless Dirac cone surface state that are 
consistent with theoretical predictions from band parity 
analysis. 



II. 7T-BONDS ON A Bi 2 Se 3 SURFACE 

We first briefly review the hybrid orbital model for 
Bi2Se3 and the existence of 7r-bonds on the (lll) r sur- 
face using the primitive rhombohedral axes (i.e., the 
(001) h in the hexagonal axes)i Considering the num- 
ber of outer-shell electrons of Bi([Xe]4f 14 5d 10 6s 2 6p 3 ) and 
Se([Ar]3d 10 4s 2 4p 4 ) and the close-packed BiSe6-SeBi6 oc- 
tahedra, chemical bonding can be described by six hy- 
bridized sp 3 d 2 orbitals to satisfy the correct coordina- 
tion numbers with two electrons per a-bond, as shown 
in Fig. [TJ Both the Bi and Se2 atoms (in the central 
Se layer of the quintuple layers) can participate in the 
six a-bonds required for octahedral shape coordination 
to create a well-defined quintuple layer. However, the 
Sel atom (the Se atom near the van der Waals gap or 
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FIG. 1. (color online) The proposed orbital hybridization for a 
quintuple Bi2Se3 layer with the actual crystal structure shown 
on the top left ( Sel (red), Bi (purple) and Se2 (orange)) and 
the Sel surface shown on the top right. The green dashed lines 
represent the six hybridized sp 3 d 2 channels with an octahedral 
shape, in which cr-bonds are formed between the Bi and Se 
atoms, but the three unpaired 7r-electrons of Sel at the van 
der Waals gap or surface may form a 7r-bond trimer with three 
of the six Sel atoms in the same layer. 

surface) has three electrons that remain unpaired after 
hybridization. To reduce Coulomb repulsion, the maxi- 
mum separation between the three unpaired electrons on 
the Sel atom corresponds to a planar configuration with 
120° angles between the electrons, i.e., a trimer, where 
three out of the six neighboring Sel atoms in the same 
layer can participate in three 7r-bonds as shown in the 
inset of Fig. [TJ We will refer to the three symmetrically 
distributed 7r-bonds as a 7r-bond trimer, and a 7r-electron 
is used to denote the unpaired electron resulting from the 
breaking of a 7r-bond. 

The 7r-bond trimer on the Sel surface of Bi2Se3 has 
two equivalent states within a six-fold symmetry coordi- 
nation, as illustrated in the inset of Fig. [TJ Because the 
7r-bond is actually formed by an extremely weak instanta- 
neous attractive force due to minimum orbital overlap, 8 
entropy is maximized by breaking and re-forming the tt- 
bonds to access the true ground state that minimizes 
the Gibbs free energy AG = AH - TAS. Clearly, the 
two equivalent states of the 7r-bond trimer in space and 
time correspond to a conjugated system of two quantum 
mechanical ensembles without a stability preference^ A 
local energy exchange mechanism between the 7r-bond 
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FIG. 2. (color online) (a) The diamond crystal structure of 
a-Sn, where all Sn atoms are in identical tetrahedral coor- 
dination, (b) The (001) surface obtained by breaking two 
of the four cr-bonds, creating two 7r-electrons per surface Sn, 
is shown. The two 7r-electrons are 180° apart and form a 
7r-bond dimer with two identical configurations, i.e., one con- 
jugated system is formed for two equivalent states, c) The 
hybrid orbital model for a-Sn is shown, where sp 3 hybridiza- 
tion is satisfied for Sn ([Kr]4d 10 5s 2 5p 2 ) with four neighboring 
Sn atoms in tetrahedral coordination. 

trimer conjugated system and six- way Se-Se connection 
on the Bi2Se3 surface (inset of Fig. [1]) has been pro- 
posed to rationalize the unique Bi2Se3 surface conduc- 
tion mechanism analogous to a dynamic three-way cross- 
bridge connecting six ports. The proposed surface con- 
duction mechanism for graphene and Bi2Se3 surfaces is 
drastically different from the traditional view of metal- 
lic conduction that is based on the diffusion of itinerant 
electrons in the conduction band. 7 



III. a - Sn: tt-BOND DIMERS IN FOUR-FOLD 
SYMMETRY COORDINATION 



Topological insulators with a four-fold symmetry in 
2D band projection have also been predicted and subse- 
quently found, such as the strained a-Sn and HgTe/CdTe 
quantum well with a zinc-blende structure ] 12 i 13 Fig. [2] 
shows the diamond structure of a-Sn. The tetrahedral 
coordination for each Sn atom ([Kr]4d 10 5s 2 5p 2 ) requires 
sp 3 hybridization using the four outer shell electrons in 
both the 5s and bp orbitals. A (001) surface cut breaks 
two of the four cr-bonds, leaving two unpaired 7r-electrons 
on the surface that can form two 7r-bonds with two of 
the four neighboring Sn atoms on the same (001) sur- 
face, as shown in Fig. [2fb). To avoid Coulomb repul- 
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FIG. 3. (color online) (a) Parity and handedness for chirality 
is shown in 2D. (b) A 7r-bond trimer in a six- fold symmetry 
coordination is shown, (c) A 7r-bond dimer in a four-fold 
symmetry coordination is shown, (d) A lone 7r-bond in a 
three- fold symmetry coordination is shown. 

sion, the maximum separation between the 7r-bonds cor- 
responds to an angle of 180° between the bonds, i.e., a 
7r-bond dimer. There are two equivalent states for each 
7r-bond dimer in space and time. Because the 7r-bond is 
extremely weak, being induced by transient attractions 
through minimal orbital overlap, it can be easily per- 
turbed by an entropy increase to lower the Gibbs free 
energy and access the true ground state. We can thus 
identify the 7r-bond dimer in four-fold symmetry coordi- 
nation on the surface as a conjugated system, in a simi- 
lar manner as for the 7r-bond trimer in six-fold symme- 
try coordination on the Bi2Se3 (lll) r surface discussed 
abovei The same cross-bridge model for surface conduc- 
tion could also be applied to the strained a-Sn as a topo- 
logical insulator when the 7r-bond dimer is semi-localized 
to maintain a dynamic conjugated system without pre- 
ferential localization or itinerancy. 

IV. PARITY SYMMETRY BREAKING OF 
7T-BONDS IN 2D 

To date two main relevant types of high symmetry sur- 
face atomic arrangements have been found in the topolog- 
ical insulators with either six- or four-fold symmetry. The 
first type is found in the (00 1)^ plane of layered Bi2Se3, 
Bii-^Sb^, and in the (111) plane of a rocksalt structure, 
such as SnTe^ii^ while the second type is found in the 
(001) plane of zinc-blende structures, such as strained 
HgTe and a-S m 12 ! 13 We can identify 2D parity symme- 
try for both a 7r-bond trimer in six-fold symmetry co- 
ordination and a 7r-bond dimer in four-fold symmetry 



coordination, as illustrated in Fig. 03 While parity trans- 
formation is related to the inversion of the spatial coor- 
dinates, the x- and y-axes cannot be inverted simultane- 
ously to obtain an inverted image in 2D, i.e., the determi- 
nant of the parity operation P must be -1, as illustrated 
in Fig. [3ja). If this condition is not satisfied, flipping 
the signs of both the x- and y-axes only produces a 180° 
rotation without mirror inversion. The parity transfor- 
mation can test for chirality by determining whether a 
mirror image is obtained following the parity operation. 
Chirality is confirmed by the breaking of parity symme- 
try. We can select between the two high-symmetry states 
by rejecting a mirror line that results in identity transfor- 
mation such that parity symmetry is broken for surfaces 
with 7r-bond-trimers in six-fold symmetry coordination 
and 7r-bond dimers in four-fold symmetry coordination, 
as illustrated by the mirror lines in Fig. G^b)-(c). In fact, 
similar parity symmetry breaking is also observed in the 
honeycomb graphene structure as shown in Fig.[3jd), for 
a properly chosen mirror linei We submit that 2D parity 
symmetry breaking is consistent with the unique charac- 
ter of surface chiral electrons with E(k) dispersion in a 
Dirac cone shape for a topological insulator. 



V. SURFACE SPIN POLARIZATION AND 
PREFERRED CHIRALITY 

The existence of chiral fermions in topological insula- 
tors in edge state (in 2D) and surface state (in 3D) with 
spins normal to the direction of their lattice momentum 
has been predicted by the quantum spin Hall effect i The 
locked spin momentum and tt Berry's phase in Bi2Se3 is 
regarded as one of the defining properties of topological 
insulators, as confirmed by spin-ARPES. 4 Chirality for a 
7r-bond trimer in six-fold symmetry coordination and a 
7r-bond dimer in four-fold symmetry coordination must 
be closely connected with the surface spin polarization. 
Generally, right-handed (RH) or left-handed (LH) chiral- 
ity is expected to be statistically equal, with no prefer- 
ence based on entropic considerations. However, the dis- 
covery of surface spin polarization strongly suggests that 
RH and LH chirality is not balanced for 7r-bond trimers 
and dimers on a topological insulator surface. 

We submit that a novel surface spin-orbit coupling re- 
sults from the chirality imbalance that follows 2D par- 
ity symmetry breaking such that the 7r-electrons effec- 
tively behave as if they were in a non-zero magnetic 
field generated by the uncanceled orbital angular mo- 
mentum, analogous to classical mechanics. This sub- 
tle effect must be a quantum phenomenon arising from 
entropic considerations for a dynamic conjugated sys- 
tem of two energetically equivalent ensembles of differ- 
ent 7r-bond arrangements. Nature clearly provides an 
even subtler mechanism for the surface 7r-electrons to ac- 
cess the true ground state through an additional mag- 
netic energy gain (—HAM) in the Gibbs free energy of 
AG = AH - TAS - HAM, in addition to the entropic 




FIG. 4. (color online) Three typical surface atomic pack- 
ings with various coordination numbers: (a) four-fold square 
packing, (b) six-fold hexagonal packing, and (c) three-fold 
honeycomb packing. 

contribution of —TAS. 

Magnetic impurities, such as Mn, have been shown 
to be capable of opening the gap in the Dirac cone for 
Bi2Se3, 4 which can also be explained intuitively using our 
proposed 7r-bond model. Mn doping effectively destroys 
2D parity symmetry by adding a strong local magnetic 
field to polarize the electron spins, creating a strong and 
random handedness preference. In other words, the sys- 
tem initially requires a very small gain in the effective 
magnetic energy —HAM to compete with the entropy 
increase —TAS from the slight chirality imbalance, how- 
ever, the strong local field created by the large Mn mag- 
netic moment breaks the subtle imbalance and shifts the 
system toward a new ground state within the correlation 
length of the magnetic moment, which does not allow for 
the formation of a dynamic 7r-bond conjugated system 
when the handedness is biased beyond being balanced 
by the entropy. 



VI. 7T-BOND MODEL PREDICTIONS FOR 
TOPOLOGICAL INSULATORS 

A general rule in crystal formation is that the num- 
ber of bonds is maximized for the greatest crystal en- 
ergy gain. Surface atomic packings typically fall into one 
of the three categories almost exclusively as a result of 
atomic size matching and outer shell electron bonding 
in 3D, as illustrated in Fig. HJ The existence of 7r-bond 
trimers in six-fold symmetry coordination and 7r-bond 
dimmers in a four-fold symmetry coordination is a com- 
monly found feature in topological insulators. The hon- 
eycomb structure that has been predicted for compounds 
such as LiAuSe could be closely related to a 7r-bond con- 
jugated system in three- fold symmetry coordination^ A 
typical signature of 7r-bonds on the surfaces of these three 



classes of topological insulators is 2D parity symmetry 
breaking, as illustrated in Fig.O In fact, we submit that 
2D parity symmetry breaking in a 7r-bond conjugated 
system is the only criteria by which a material may be 
considered to be topologically trivial or nontrivial. 

The necessary condition for the existence of a dynamic 
7r-bond conjugated system in topological insulators is the 
requirement of weak covalent bonds in the material. For 
example, nearly all of the verified topological insulators 
fall into the category of narrow band gap semiconduc- 
tors, the main character of which is a low melting point 
resulting from weak covalent bonding. If the material is a 
semimetal with a small band overlap, the band gap can be 
adjusted through proper atom substitution, such as tun- 
ing Bi into a narrow band gap semiconductor so that the 
semiconductor acts as topological insulator Bii-^Sb^, 9 
In contrast, for insulators with covalent bonds that have 
a slight polarity to ionic character due to a relatively large 
EN difference, atomic substitution can narrow the band 
gap by EN adjustment; for example, the EN difference 
of PbTe with a rocksalt structure (EN of Pb=2.33 and 
EN of Te=2.1) may be tuned with Sn substitution (EN 
of Sn=1.96). 16 We conclude that the role of the strain in 
tuning a topological state is actually a side effect of the 
atomic substitution because the strain is introduced by 
the unavoidable mismatch in atomic sizes through sub- 
stitution, although the goal is electronic band gap size 
tuning. 

Based on our proposed 7r-bond model, HgTe with a 
zinc-blende structure is likely to be topologically non- 
trivial due to a low melting point from weak covalent 
bonding (EN of Hg=2.0 and Te=2.1), and the existence 
of surface 7r-bond dimers in four-fold symmetry coordi- 
nation on (001) planes, in agreement with the prediction 
by Chadov et al^ However, we expect AmN with a rock- 
salt structure to be topologically trivial^ mainly due to 
the particularly large EN difference between the Am and 
N atoms. The existence of the dynamic 7r-bond conju- 
gated system on topological insulator surfaces requires 
a low surface energy to accommodate the entropic fac- 
tor on lowering the Gibbs free energy. The energies of 
the itinerant electrons in metals and the unpaired elec- 
trons resulting from the strong covalent bonds in a wide 
band gap insulator being broken are clearly too high to 
maintain even 7r-bonds on the surface. 



VII. CONCLUSIONS 

In summary, we have applied a hybrid orbital model to 
typical topological insulators for Bi2Se3 and a-Sn based 
on crystal symmetry, surface atomic coordination, and 
the outer-shell electrons. The identification of 2D parity 
symmetry breaking leads to a reasonable explanation for 
the chiral character of surface electrons, which is implied 
by the quantum spin Hall effect, and this chirality im- 
balanced can be used to interpret the spin polarization 
found in the surface band of a Dirac cone. It is clear 
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that the 7r-bond is key to understanding the surface elec- 
tronics, as it has been repeatedly found in the ideal 2D 
graphene electron gas and topological insulator surfaces 
with similar Dirac cones. 
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